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Abstract 

Effective bending moduli of a fluid membrane are investigated by means of the transfer-matrix 
method developed in our preceding paper. This method allows us to survey various statistical 
measures for the partition sum. The role of the statistical measures is arousing much attention, 
since Pinnow and Helfrich claimed that under a suitable statistical measure, that is, the local 
mean curvature, the fluid membranes are stiffened, rather than softened, by thermal undulations. 
In this paper, we propose an efficient method to observe the effective bending moduli directly: 
We subjected a fluid membrane to a curved reference plane, and from the free-energy cost due 
to the reference-plane deformations, we read off the effective bending moduli. Accepting the 
mean-curvature measure, we found that the effective bending rigidity gains even in the case of 
very flexible membrane (small bare rigidity); it has been rather controversial that for such non- 
perturbative regime, the analytical prediction does apply. We also incorporate the Gaussian- 
curvature modulus, and calculated its effective rigidity. Thereby, we found that the effective 
Gaussian-curvature modulus stays almost scale- invariant. All these features are contrasted with 
the results under the normal-displacement measure. 
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I. INTRODUCTION 



Amphiphilic molecules in water segregate spontaneously into flexible extended surfaces 
called fluid (lipid) membranes IlL [M . The fluid membranes are free from both surface tension 

n n 

and shear modulus, and the elasticity is governed only by bending rigidity [3|, |4J. The 
Hamiltonian is given by the following form, 

H = JdA (~ J 2 + kK\ . (1) 

The mean curvature J is given by the summation of two principal curvatures J = c\ + C2, 
whereas the Gaussian curvature K is given by their product K = c\Ci. The corresponding 
two moduli k and R are called bending rigidity and Gaussian-curvature modulus, respec- 
tively. The integration J dA extends over the whole membrane surface. The Gaussian- 
curvature term governs the global structure of the membranes, because the term / dAK 
measures a topological index; for instance, for a vesicle with rih handles, such an identity 
/ dAK = 47r(l — rth) holds (Gauss-Bonnet theorem). 

In spite of its seemingly simple expression, it is very hard to treat the Hamiltonian ((TJ 
by analytical methods. As a matter of fact, when written in terms of an explicit parame- 
terization, the Hamiltonian becomes ugly; see Eqs. (??)-(??) mentioned afterwards. Hence, 
owing to the thermal undulations and the mutual interactions, it is expected that the effec- 
tive bending moduli are modified effectively for macroscopic length scales. In orde_r to_ clarify 
this issue, numerous renormalization-group analyses have been reported so far 
For the bending rigidity, the following renormalization-group equation has been obtained; 

7 rji 

with renormalized bending rigidity k', temperature T, and the number of decimated 
molecules M. Literature agrees that the numerical prefactor in the above equation is a = 3. 
[A more detailed account of the historical overview would be found in Ref. [jj.] Because of 
a > 0, the effective bending rigidity is reduced by thermally activated undulations. This 
conclusion might be convincing, because the membrane shape itself should be disturbed by 
the thermal undulations. As a matter of fact, it has been known that the orientational 
correlation is lost for long distances 1(3]. It is quite natural to anticipate that membranes 
become flexible for length scales exceeding this correlation length. 
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Recently, however, Pinnow and Helfrich 0, obtained a remarkable conclusion a 
— 1(< 0). The key ingredient of their new argument is that they considered the role of 
measure factors for the partition sum. They insist that the local mean curvature J should 
be the right statistical measure rather than other measures such as the normal displacement 
h and the local tilt angle 9. [The normal displacement h has been used as a standard measure. 
We will explain the /i-based parameterization afterwards.] After an elaborated calculation of 
the variable replacement h — > J and succeeding renormalization-group analysis, the authors 
reach the conclusion of a = —1. Moreover, as for the Gaussian-curvature modulus, they 
insist that the effective modulus k' should remain scale-invariant; 

R' = R. (3) 

This conclusion again contradicts the common belief that R would be enhanced for macro- 
s co pl e length S ea!e S ; llamdy , k' — K + ( W ) ,nM gQ. This enh_ t sign* the 
topological instabilities. 

The developments mentioned above all stem from the approximation that the membrane 
is almost flat, and the thermally excited undulations are extremely small. In our preceding 



paper [12|, we developed an ab initio simulation scheme in order to study the thermodynam- 
ics of a fluid membrane beyond such perturbative level. As a demonstration, we calculated 
the transformation coefficient dn'/ dn, which yields the direction of the renormalization-group 
flow through coarse-graining. We found that the renormalization-group flow is influenced 
significantly by the choice of the statistical measures. In fact, under the mean-curvature 
measure, we observed an indication that the effective bending rigidity flows toward the 
large-K direction. We did not include the R term in the preceding work. 

In this paper, extending the preliminary analysis, we propose an efficient method to 
observe the effective bending moduli directly: We subject a fluid membrane to a curved 
reference plane. From the free-energy cost due to the reference-plane deformations, we read 
off the effective bending moduli. Our data indicate definitely that under the mean-curvature 
measure, the membrane stiffening occurs even for the case of very flexible membrane k < 
1. This result supports the picture that the stiffening, contrary to our naive expectation, 
is driven by the thermal fluctuations. We also calculate the effective Gaussian-curvature 
modulus, and show that it remains almost scale- invariant. Again, the result is in good 
agreement with the Pinnow-Helfrich claim. 
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It has to be mentioned that the Monte-Carlo method has been utilized successfully in 
the studies of membranes and vesicles |2], For the Monte-Carlo method, however, a teth- 
ered (polymerized) membrane [13] rather than a fluid membrane is more suited, because a 
membrane is implemented in a computer as an assembly of molecules and junctions bear- 
ing close resemblance to a tethered membrane. (Note that because of the absence of shear 
modulus, fluid membranes should have no internal structure.) However, Gompper and Kroll 
succeeded in simulating fluid membranes by the Monte-Carlo method, allowing reconstruc- 



tions of junctions during the simulation 



14j . They observed the topological instabilities 



with respect to the variation of temperature and membrane concentration. In fairness, it 
has to be mentioned that their Monte-Carlo data indicate softening for lipid vesicles. 

The rest of this paper is organized as follows. In the next section, we describe our 



i. We also explicate 
. In Sec. ??, we 



new approach to the direct observation of the effective bending modu 
the outlines of the numerical method which was reported previously 
present the numerical results. We are mainly concerned in the case of the mean-curvature 
measure. For a comparison, we will also present the data calculated under the normal- 
displacement measure. With the first-principle simulation method, we will show that the 
scenario advocated by Pinnow and Helfrich holds true even for the non-perturbative regime. 
In the last section, we give summary and discussions. 

II. LEGENDRE TRANSFORMATION AND THE EFFECTIVE BENDING MOD- 
ULI 

In this section, we explain the idea for calculating the effective bending moduli. Math- 
ematical formulas necessary in the succeeding numerical simulations are derived. We start 



12 1 . It is im- 



with recalling the outlines of the transfer-matrix method proposed previously 
portant to recognize the outlines of the transfer-matrix construction, because it elucidates 
the underlying physical idea of our effective-bending-moduli calculation. In short, it has to 
be recognized that our membrane should be classified into the "open framed membrane" in 
the category of Ref. J^. 



4 



t of the microscopic constituents realizing the Hamiltonian Hence, 
1^ | . an alternative, in a sense, rather simple-minded, approach to the 



A. Transfer-matrix formalism: A brief remainder of Ref, 

As mentioned in Introduction, the fluid membranes are free from shear modulus. This 
fact tells that the fluid membranes should have no internal structure. Hence, it is by no 
means fruitful to thin 
we proposed in Ref. 

fluid membrane: We constructed the transfer matrix directly from the Hamiltonian {IJ). In 
the construction, we managed several discretizations which we explain below. 

As noted in Introduction, the membrane shape is parameterized by the normal (trans- 
verse) displacement h(x, y) from a base (reference) plane. The variables x and y denote the 
Cartesian coordinates on the reference plane. In term of this displacement field h(x,y), the 
mean curvature and the Gaussian curvature are parameterized explicitly as follows |15| : 

_ {&jh + d 2 y h)(l + {d X hf + {dyhf) ~ 2d X hdyhd X dyh ~ ^h(d x h) 2 ~ d'hjdyh) 2 
[X,V) (1 + {d X hf + {dyhff/ 2 

(4) 

and, 

respectively. Similarly, the infinitesimal area dA is given by, 

dA={l + (d x hf + {d y h) 2 ) 1/2 dxdy. (6) 

Putting them together into the Hamiltonian (0), we arrive at the explicit representation 
in terms of the displacement field h(x,y). Now, we are led to a two-dimensional scalar- 
field theory with considerably complicated interactions. It is notable that the theory has 
even no obvious perturbation parameter. This fact motivated us to develop a first-principle 
simulation technique. 

We put the theory on a square lattice with lattice constant a; see Fig. ?? (a). Accordingly, 
the field variables are now indexed by integer indices; namely, h(x,y) — > hij. Hereafter, we 
set the lattice constant as the unit of length; namely, a = 1. Our theory has the translational 
invariance of h — > h + Ah, and thus, the absolute value of h is meaningless. Therefore, it 
is sensible to use the link variable s = adh rather than h; note that s denotes the step 
associated at each link; see Fig. ??. We are led to the dual lattice model. Thereby, for each 
shaded plaquette of this dual lattice (originally a vertex spanned by four adjacent links), 
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the following local statistical weight is associated; 



p(si, s 2 , s 3 , Si) = exp(-H), 



(7) 



with, 



Ti = dA(s x , s 2 , s 3 , s 4 )(^ J (si, s 2 , s 3 , s 4 ) 2 + K,K(s h s 2 , s 3 , s 4 )), 



(8) 



where J(s\, s 2 , s 3 , s 4 ), X(si, s 2 , s 3 , s 4 ) and cL4(si, s 2 , s 3 , s 4 ) are to be replaced |l2j with the 
finite-difference versions of the differential forms (??)-(??). See Fig. ?? (b) for the definitions 
of {s a }. We have set ksT = 1, because this factor can be absorbed into the bending moduli 
k and R. Moreover, we should introduce yet another "statistical weight" for each open 
plaquette so as to impose the constraint of rots = 0; the gradient field should be rotationless. 
That is, 

A(sj, s 2 , s 3 , s 4 ) = <5 Sl+S2 _ S3 _ S4j0 - (9) 

To summarize, we are led to the dual lattice model with the step variable s. There 
are two types of statistical weights p{s\, s 2 , s 3 , s 4 ) (??) and A(si, s 2 , s 3 , s 4 ) (??), which are 
arranged in the checkerboard pattern. Likewise the transfer matrix is constructed as a strip- 
like segment shown in Fig. ?? (c). It is a good position to mention a number of remarks: 
First, the above theory takes the displacement variable hij as the statistical measure. As 
noted in Introduction, the mean-curvature statistical measure is considered to be physically 
sensible. The conversion of the statistical measure is achieved by the redefinition of the 
statistical weight; namely, 



p(si, s 2 , s 3 , s 4 ) — > p(si, s 2 , s 3 , s 4 ] 



\ 



n 

a=l 



dJ(s 1 ,s 2 ,s 3 ,s 4 ) 



(10) 



The square root is intended to take the geometrical mean, because each step variable s a 
is sheared by an adjacent plaquette as well. Second, the step variable is discretized as 
Si = S s (i — N s /2 — 0.5) with i = 1, • • • , iV s . The unit of step S s is determined self-consistently 
during the simulation by 5 S = RJ {sf). This step-variable discretization is an influential 
factor concerning the reliability of the present simulation, and its performance was checked 
previously . Third, in order to diagonalize the transfer matrix, we utilized the the 

nnnn 
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density-matrix renormalization group 



so as to investigate the highly-correlated systems such as the Hubbard models and the spin 
chains. Later on, it was extended to the field of soft materials such as the lattice vibrations 



2jJ, |22j, |23| , the quantum string j^, 25], and the bosonic systems j^El- We repeat density- 
matrix renormalization one after another so as to reach sufficiently long transfer-matrix strip 
length . The number of states retained for a renormalized block is an important technical 
parameter which is denoted by m. 

One may wonder that the strong deformations increase the area of the piece of membrane 
considered to such an extent that the emerging additional degrees of freedom for additional 
molecules cannot be ignored. According to the idea of Helfrich, however, the correct sta- 
tistical measure is the local mean curvature that has noting to do with the "molecules." 
Therefore, our treatment is justified even for such strong deformations as long as we accept 
his idea. Strictly speaking, we are not dealing with such "molecules" at all; rather, we had 
just discretized the real space so as to form a square network, and the vertices are not to be 
regarded as molecules. In fact, the fluid membrane should have neither internal structure 
nor fixed connectivity, and it is not fruitful to think of microscopic constituents, and regard 
them as the degrees of freedom. 

Finally, in the above, we have postulated the presence of a reference plane from which all 
undulations are created, and we imposed no restriction to the total area of the membrane. 
Hence, in terms of the category of Ref. [2j, the membrane is to be classified into "open 
framed membrane;" see Ref. 29] as well. The free energy per unit area of the reference 
plane is the naturally observable quantity. This leads to the idea that the effective bending 
moduli would be measured from the increase of the stress energy due to the reference-plane 
deformations. We will pursue this idea in the next subsection. 



B. Legendre transformation and the effective bending moduli n e jj and reifV 

In this subsection, we explain our approach to the effective bending moduli. There are 
two types of bending moduli such as k and k; see Eq. First, we explain the way to 

calculate the effective bending rigidity. We introduce the following Hamiltonian with an 
additional term; 

H C = H - £ CJ, = £(^eL4, + RKidAi - CJi). (11) 

The index % runs over the shaded plaquettes of Fig. ??; The quantities Jj, K i: and dAi are 
the same as those in Eq. (??), but possessing the plaquette index i now. Note that the 
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additional term is not lumped together with the factor cL4j, because our aim is to calculate 
the stress energy with respect to the reference plane rather than the membrane surface itself. 
The additional term breaks the symmetry of the mean curvature J <-> — J linearly. In other 
words, the membrane is forced to bend so as to possess a non-vanishing spontaneous mean 
curvature. Hence, with respect to the stress energy due to the reference-plane deformation, 
we are able to observe the effective bending rigidity. That is the basic idea of our approach. 
Such idea was featured in an analytical treatment as well [jj. By the way, our aim is to put 
forward this idea to an actual computer simulation. 

As is well-known, the above idea is best formulated by the Legendre transformation; 

G(j) = F(C) + Cj, (12) 

with dF/dC = —j. Here, F(C) denotes the free energy of the Hamiltonian (??) per unit 
cell (one shaded plaquette). G(j) is the desired Legendre-transformed free energy, which 
is a function of the spontaneous mean curvature j. Our concern is to obtain the effective 
bending rigidity, 



K eff 



d 2 G 



With use of the well-known identity (d 2 G/dj 2 )(d 2 F/dC 2 ) = —1, we obtain the expression 
for the effective bending rigidity; 

= -V^- (14) 

Let us address a number of remarks: First, the free energy F is readily accessible by the 
transfer-matrix calculation. Hence, the above formula is suited to our computer simulation. 
The remaining task is the numerical differentiation. We had adopted "Richardson's deferred 
approach to the limit" algorithm explicated in the textbook |2j|. Second, we started from 
the Hamiltonian (??), which is defined on the reference plane rather than the original fluctu- 
ating membrane surface. Therefore, the effective bending rigidity corresponds to the elastic 
modulus with respect to the reference-plane deformation. That is precisely what we sought. 

The expression of Eq. (??) may seem to be exceedingly formal. In fact, it may be 
unclear how the interaction of the curvature with the undulations is taken into account. 
It is noteworthy that the Hamiltonian contains the symmetry breaking term CJi and we 
have to evaluate the free energy in the presence of it. Just like the two-dimensional Ising 
model with the external field that has not yet been solved exactly, such problem with the 
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symmetry breaking term is far from being trivial in itself, and generally, it contains valuable 
informations such as the interaction between the background curvature and the thermal 
undulations. In fact, Eq. (??) states that the second derivative with respect to C yields 
the effective rigidity. That is the underlying idea behind the formal expression of Eq. (??). 
We stress that such an additional term is readily tractable by our simulation method, and 
owing to this advantage, we are able to access the effective rigidity in a quite straightforward 
manner. 

We will introduce another effective bending rigidity: Through the coarse-graining de- 
picted in Fig . ??, we obtain a coarse-grained lattice and the corresponding smeared cur- 



vature J 12|. after rescaling the unit of length \/2a — > a and reexpressing the formulas in 
terms of the rescaled (smeared) quantities [we put "for them], we obtain the effective rigidity 
for the coarse-grained membrane; 



Note that two unit cells are renormalized into one coarse-grained unit cell. In other words, 
two molecules are renormalized into one decimated molecule, and hence, the parameter M 
in Eq. © should be M = 2. 

Let us turn to the Gaussian-curvature modulus. In this case, we incorporate the following 
additional term coupling to the Gaussian curvature linearly; 

H D = H-J2DK t = Y,{^dA t + RKMi ~ DK { ). (16) 

i i 1 

Similar to the above, this leads to the following Legendre transformation, 

G(k) = F(D) + Dk, (17) 

with dF/dD = —k. Our concern is to obtain the effective Gaussian-curvature modulus, 
which, in other worlds, the effective symmetry breaking term with respect to the K field. 
Postulating that the effective free energy G(k) is a quadratic polynomial in terms of k, we 
obtain the following expression, 

- JJL/^L (1 8) 

Again, the similar idea applies to the coarse-grained lattice. Hence, we obtain, 



K 



(L) _ OF ,d 2 F 
eff ' 3D' dD 2 
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k/^r- ( 19 ) 



We complete preparing the mathematical formulas for the effective bending moduli. We 
apply these formulas to the computer simulation in the next section. 



III. NUMERICAL CALCULATION OF THE EFFECTIVE BENDING MODULI: 
ROLE OF THE STATISTICAL MEASURES 



K 



In this section, we explore the effective bending moduli [k^'P, Eqs. (??) and (??); and 

(S,L) 
eff ' 



Eqs. (??) and (??)] of a fluid membrane by means of the transfer-matrix method 



explained in Sec. ??. The transfer matrix is diagonalized 12] by means of the density- 



matrix renormalization group 



17 



18. 
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. Making a comparison between the results under 



the mean-curvature and the normal-displacement measures, we will elucidate the role of 
statistical measures. As mentioned in the above section, we have fixed the temperature 
(ksT = 1), because this factor can be absorbed into the redefinition of the bending moduli. 

Here, we shall outline some technical points, that are relevant to the simulation precision. 
[Detailed account of the simulation algorithm is presented in Ref. [l^.] We repeated forty 
renormalizations for obtaining each plot; namely, the strip length of the transfer matrix 
extends to L — 80. The technical parameters m, N s and R are indicated in each figure 
caption; see Sec. ?? for the meanings of these technical parameters. These parameter values 



are equivalent or even improved to those used previously 
the simulation is maintained. 



121 ] . Therefore, the reliability of 
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A. Effective bending rigidity 

In this subsection, we focus our attention on the effective bending rigidity. For that 
purpose, for the time being, we drop the Gaussian-curvature term, which is considered in 
the next subsection; namely, we set R = 0. In Fig. ??, we plotted the effective bending 
rigidity dftp for various bare rigidity k. The moduli K^fX and denote the effective 
bending rigidities for smaller and longer length scales, respectively; see Eqs. (??) and (??). 
Here, we have accepted the local mean curvature J as for the statistical measure; recent 
through discussion |<J, [ll| insists that this statistical measure should be the right one. From 



the plot, we see that the rigidities and deviate from each other as for small k. In 
the small-K regime, the membrane becomes very flexible, and so the thermal fluctuations 
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iffening 

m 



should be enhanced significantly. Hence, we see that the correction to the effective rigidity 
is actually induced by the thermal fluctuations. Moreover, we notice Krjff > «e//- Hence, 
the membrane acquires stiffness for longer length scales; namely, the membrane stiffe ning 
sets in. Such membrane stiffening was first predicted by the analytical arguments 
However, for the small-K regime, the analytical arguments are not fully justified, because 
the arguements stem from the "nearly flat approximation." On the other hand, our first- 
principle simulation does not rely on any perturbative treatment. In that sense, our data 
demonstrate very definitely that the membrane stiffening occurs withstanding the thermal 
disturbances. 

Note that the effective rigidity K K e ^ is by no means identical to the "bare" coupling 
constant k. The former is the bending elasticity with respect to the reference-plane defor- 
mations, whereas the latter is the elastic constant of the membrane surface itself. Therefore, 

they need not coincide. However, for sufficiently large k, as is seen from Fig. ??, the 

(s) 

curve tends to be parallel to the slope of the line K e f± = k, indicating that they coincide 
asymptotically for large k. 

We shall argue the relationship between the above result and our previous report jl^ . 
The regime k < 0.4, where we found a notable deviation of nf^ and k^X in Fig. ??, 
coincides with the area of the prominent dn' /dn enhancement reported in Fig. 7 of Ref. 



12j |. [Although dn' /dn does not yield direct assessment of the effective rigidity, we concluded 



that the dn' /dn enhancement should reflect the membrane stiffening.] Hence, in retrospect, 
our preceding analysis appears to capture the precursor of the membrane stiffening fairly 
correctly. 

In Fig. ??, at k ~ 0.4, there appears a singularity: The moduli vifX and k^X approach 
to each other, and for k > 0.4, they split off again. This singularity may indicate an onset of 
a phase transition. For k > 0.4, because of K^fj > ftfX, a membrane stiffening should occur. 
As a matter of fact, because of the discretization of the step variables [see Sec. ??], it is 
likely that the membrane is trapped by the most stable configuration (flat surface) for large 
k. [It is expected that the membrane becomes flat just like the smooth phase in the solid- 
on-solid model with large surface tension^ Moreover, it is well known that the correlation 
length diverges exponentially for large k [lOj]. Such long correlation length would exceed the 
capability of the numerical simulation. The large- k behavior appearing in k > 0.4 is thus 
an artifact of the numerical simulation. The membrane stiffening for large k is not intrinsic, 



11 



and is rather driven by the mechanism different from that of the small- k side. 

In Fig. ?? (a), keeping such drawback in mind, we have drawn an anticipated behavior 
of the effective bending rigidity for a wide range of k. As mentioned above, the analytical 
treatment is justified for large k. On the other hand, our first-principle simulation works 
efficiently in the other side (non-perturbative regime). For sufficiently large k, the analytical 
argument predicts the correction to the effective k such lnM/8ir = 0.027. . .. 

The correction appears to be exceedingly small to be resolved by the numerical simulation; 
see Fig. ?? as well. On the other hand, for the small-K regime, our simulation data indicate 
that the correction to the effective k increases very significantly. The amount of correction is 
comparable to the thermal-fluctuation energy; note that we have chosen fcgT as the unit of 
energy (/cgT = 1). Hence, it is suggested that for macroscopic length scales, the membrane 
acquires a considerable amount of effective stiffness, and it would look almost flat irrespective 
of the thermal disturbances. 

Let us turn to the normal-displacement statistical measure. In Fig. ??, we plotted the 

(S V) 

effective bending rigidity K e fj under this statistical measure. Notably enough, the behavior 
is quite contrastive with that of the local-curvature measure mentioned above: The larger- 
scale effective rigidity k^X is suppressed by the thermal undulations. Hence, it is shown 
that the membrane is softened effectively for macroscopic length scales. This result may 
meet our intuition, and has beenpredicted by numerous analytical arguments based on the 
normal-displacement measure . We stress that our first-principle simulation covers 

various statistical measures in a unified way. Our simulation clarifies fairly definitely that 
the choice of measure factors is vital for the thermodynamics of the fluid membrane. Again, 
for large k ~ 0.8, a signature of the membrane stiffening k^X > k^X comes up, and the 
membrane should undergo the flat phase. This behavior is a drawback of our simulation 
as mentioned above. Keeping this in mind, we have drawn a schematic behavior of the 
effective bending rigidity in Fig. ?? (b). For large k, the analytical argument predicts the 
renormalization correction k' — k — —3lnM/8n = —0.082 . . ., which is beyond the resolution 
of the present numerical simulation. As for the small-K regime, the correction is enhanced. 
However, the enhancement is not so prominent as in the case of the mean-curvature measure. 
It is almost comparable to the analytical prediction; see Fig. ?? as well. Hence, our first- 
principle simulation indicates that the analytical formula (J2J) is more or less applicable even 
for the case of the non-perturbative regime under the normal-displacement measure. 
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B. Effective Gaussian-curvature modulus dftf 

In the above, we have studied the thermal-fluctuation-induced corrections to the bending 
rigidity k. In this subsection, we incorporate the Gaussian-curvature modulus R, and look 
into its effective strength R^f ■ As n °ted in Introduction, the Gaussian-curvature- modulus 
term is related to the topological index, and hence, it governs the global structure of the 
membranes. Roughly speaking, for R > 0, the plumber 's-nightmare phase (lamellar with 
tunnel-like defects) is stabilized, whereas for R < 0, the formation of vesicles (droplets) is 
favored. In that sense, the quantity R e ^ reflects the tendencies toward the topological 
instabilities. 

In Figs. ?? and ??, we presented the effective Gaussian-curvature modulus for various 
bare R under the statistical measures of the mean curvature and the normal displacement, 
respectively. Here, k is fixed to be k = 2/ v / 2(= 0.35 . . .) and 0.4 for respective figures. First, 
let us argue the latter. [Because this case provides a prototypical example, we will argue 
it prior to the mean-curvature case.] The latter case has been studied extensively so far 
with analytical approaches 0,0]. However, because we are supposing that the membrane 
is framed by the reference plane, there emerge some characteristic features. For R « 0, we 

— (S) 

see that a large amount of effective \K e ff\ appears. This result indicates that the membrane 
undulations are dominated by the dimple-like deformations, and possibly, the membrane 
tends to form droplets. This feature is in accordant with the previous claim j^] that 
for small membrane concentration, the membranes are thermodynamically unstable to the 
dissolution into the solvents (sponge phase). On the contrary, the effective modulus for the 
longer length scale, namely, |«e//|j exhibits considerable suppression. Notably enough, it 
becomes even positive for R > 0.1: For macroscopic length scales, the membrane recovers 
its stability around R « 0, although microscopic undulations are in favor of droplets. Such 
crossover behavior is convincing, because the membrane is framed by the reference plane, 
and macroscopically, lamellar-type structure should be retained. 

As the bare modulus R decreases, the effective moduli of different length scales coincide 
at R ~ —0.6. For the region exceeding this point R < —0.6, in turn, R^ dominates R^X 
eventually. Hence, in this region, the droplet formation is favored for macroscopic length 
scales. The location of this transition point is reminiscent of that advocated by the analytical 
arguments (sii l^ . which predicts the transition point R c ~ — 10k/9 = —0.44. . .. In our 
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simulation, the planar-type morphology is assumed a priori. Therefore, in such regime 
R < k c , our numerical simulation does not cover such sponge phase. In fact, it suffers from 
pathologies such as the diverging mean deviation of the step variables. In the same way, 
for exceedingly large R > 0.1, the membrane becomes unstable owing to the topological 
instability toward the plumber's-nightmare phase. In the regions depicted in Figs. ?? and 
??, the planar-type morphology is retained, and thus our simulations are reliable. 

Second, let us turn to the case of the mean-curvature measure: In Fig. ??, we plot- 
ted i^fff ^ for various bare Gaussian-curvature modulus k. There appear some behaviors 
characteristic of the mean curvature measure: We notice that and R^X almost overlap 
each other around R ~ 0. That is, the Gaussian-curvature modulus stays almost scale- 
invariant through coarse-graining. This result supports the aforementioned analytical pre- 
diction H 11| of Eq. (JSJ). In addition to this, it is to be noted that ftfff'P exhibits a large 

_(S L) — — 

negative residual value k^A —1.4 around k ~ 0. That is, although the bare coupling n is 
turned off, the membrane undulations are dominated by the dimple-like deformations. This 
result validates the postulation by Helfrich insisting that the elementary excitations of 
the thermal undulations should be the "hat excitations" rather than the ordinary sinusoidal 
ones. Because the hat-excitation picture is the starting point of his arguments, the whole 
theoretical theoretical structure appears to be selfconsistent from our first-principle data. 



IV. SUMMARY AND DISCUSSIONS 

We have investigated the effective bending moduli, K K e Jf [Eqs. (??) and (??)] and R e jj 
[Eqs. (??) and (??)], with an emphasis on the role of the statistical measures for the 
mrtition sum. We employed the transfer-matrix method developed in our preceding paper 
1^ | . where we had reported a preliminary analysis on the effective bending rigidity via 
the transformation coefficient dn'/dn. In the present paper, we proposed the scheme to 
determine the effective bending moduli directly: We calculated the free-energy cost due 
to the reference-plane deformations, from which we read off the effective bending moduli. 
This idea is formulated in terms of the Legendre transformation, and the mathematical 
formalism is developed in Sec. ??. Based on the formulas, we carried out extensive computer 
simulations in Sec. ??. As a result, we found a clear evidence of the membrane stiffening in 
the case of the mean-curvature measure; see Figs. ?? and ?? (a). The membrane stiffening 
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was first predicted by the analytical approaches (3, , which are validated for sufficiently 
large k. Our first-principle data show that the membrane stiffening occurs even for the non- 
perturbative (small k) regime withstanding the thermal disturbances. Surprisingly enough, 
the enhancement of the effective bending rigidity copes with the thermal-fluctuation energy 
~ ksT, suggesting that the membrane would stay almost flat for macroscopic length scales. 

On the contrary, under the normal-displacement statistical measure, we found a clear 
indication of the membrane softening; see Figs. ?? (b) and ??. This fact indicates that 
the choice of measure factors is indeed significant. The correction to the effective bending 
rigidity appears to be moderate compared with that of the mean-curvature measure. In 
fact, it is almost comparable to the prediction by the analytical treatment of Eq. (J2J) even 
for small k. 

For exceedingly large rigidity, the membrane fluctuations freeze because of the exponen- 
tially diverging correlation length and the pinning potential due to the step-variable 
discretization; the membrane undergoes the flat phase eventually just like the solid-on-solid 
model with large surface tension. The appearance of such phase is a drawback of the nu- 
merical simulation, and in this respect, the simulation and the analytical treatment are both 
complementary. 

Furthermore, we incorporated the Gaussian-curvature modulus, and studied its effective 
strength. Accepting the mean-curvature measure, we found that the effective Gaussian- 
curvature moduli for different length scales overlap each other around R ~ 0; see Fig. ??. In 
other words, the Gaussian-curvature modulus stays almost scale-invariant through coarse- 
graining. This fact is in good agreement with the analytical prediction of Eq. In 
addition to this, the effective Gaussian-curvature modulus exhibits a large negative residual 
value even for zero bare modulus. This fact indicates that the membrane fluctuations are 
governed by the dimple-like deformations that should be scale-free. This observation again 
supports the scenario of Ref. j^f insisting that the thermal undulations should be decomposed 
into the hat excitations. 

rave investigated 



111 ]. Thereby, we 



To conclude, by means of the first-principle simulation technique, we 
the series of analytical predictions advocated by Pinnow and Helfrich 0, 
found that these predictions hold true even for the non-perturbative regime. In particular, 
the hat-excitation picture, which is the very starting point of their argument, is validated 

— (S L) 

by our R e t'f data. Therefore, the postulation and its deductive hypotheses turn out to be 
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fairly consistent as a whole. Hence, it is very likely that the mean curvature is indeed a 
physically sensible statistical measure for the partition sum. 
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(c) 

FIG. 1: (a) On the square lattice, we consider scalar field hij denoting normal displacement of a 
membrane with respect to a reference plane. Step variable (gradient field) s = adh is defined at 
each link, (b) The local statistical weights p [Eq. (??)] and A [Eq. (??)] are represented by shaded 
and open squares, respectively. The statistical weight p has a variant so as to take account of other 
integration measure such as the local mean curvature (??). (c) From these local statistical weights, 
we construct a strip whose row-to-row statistical weight yiel 
transfer matrix is diagonalized with the DMRG method 



the transfer-matrix element. This 
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FIG. 2: Real-space decimation procedure. From the decimation, coarse-grained curvatures J and 



K are constructed 



moduli; namely, k^X (??) and k^X (?? 



12|. J and -K" are used so as to obtain the corresponding effective bending 

(L) 
eff 
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FIG. 3: Effective bending rigidity K e ^ is plotted for various bare bending rigidity k and the 
fixed Gaussian-curvature rigidity R = 0. We have accepted the local curvature as for the statistical 
measure. The simulation parameters for each symbol are (+) m = 15, N s = 8, and R = 0.9; (x) 
m = 15, N s = 8, and R = 0.8; and (*) m = 15, N s = 9, and R = 0.7. Because of t^X > i^fff, we 
see that the membrane is stiffened effectively for macroscopic length scales. The stiffening in the 
large-K side, n > 0.4, may be the artifact of the numerical simulation due to the pinning potential 
of discretized step variables; there emerges the smooth phase just like the solid-on-solid model with 
large surface tension. 
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FIG. 4: Schematic drawings of reifV anticipated from our first-principle data [Figs. ?? and ??], 
and the analytical result justified for sufficiently large n [Eq. ©]. (a) The mean curvature is 
accepted as for the statistical measure, (b) The normal displacement is accepted for the statistical 
measure. 
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FIG. 5: Effective bending rigidity n^ff is plotted for various bare bending rigidity k and the 
fixed Gaussian-curvature rigidity R = 0. We have accepted the normal displacement as for the 
statistical measure. The simulation parameters for each symbol are (+) m = 13, N s = 9, and 
R = 0.55; (x) m = 10, iV s = 11, and R = 0.45; and (*) m = 11, N s = 10, and R = 0.5. Because 
of ft^j < ^ejj, we see that the membrane is softened effectively for macroscopic length scales. 
The stiffening in the large-K side, k > 0.8, may be the artifact of the numerical simulation due to 
the pinning potential of discretized step variables; there emerges the smooth phase just like the 
solid-on-solid model with large surface tension. 
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FIG. 6: Effective Gaussian-curvature modulus R e f't [Eqs. (??) and (??)] is plotted for various 
bare modulus R and the fixed bending rigidity k = 2/y/2(= 0.35...). We have accepted the 
local curvature as for the statistical measure. The simulation parameters for each symbol are (+) 
m = 15, N s = 7, and R = 1; and (x) m = 15, N s = 8, and R = 0.9. We see that the effective 
modulus stays almost scale-invariant around ~ ^ through coarse-graining, confirming the 

validity of the analytical prediction [Eq. Q]. Moreover, we notice that K^f'P exhibits a large 
negative residual value even for zero bare R = 0. This fact reflects that the membrane undulations 
are dominated by the formation of "hat excitations" |9J. 
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FIG. 7: Effective Gaussian-curvature modulus R e A [Eqs. (??) and (??)] is plotted for various 
bare modulus R and the fixed bending rigidity n = 0.4. We have accepted the normal displacement 
as for the statistical measure. The simulation parameters for each symbol are (+) m = 14, N s = 8, 



and R = 0.6: and fx) m 



11, and R = 0.45. We see that \k 



(L) 



v e jj\ is suppressed around 
R ~ 0. This fact tells that for macroscopic length scales (because of the restriction of the reference 
plane) the planar-type morphology is favored; namely, R is irrelevant in the infrared limit. For 
R < —0.6, in turn, |/Cg«| is enhanced eventually, suggesting that the membrane tends to dissolve 
into the solvent (droplet phase). 
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